Many-body effects in van der Waals-Casimir interaction between graphene layers 
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Van der Waals-Casimir dispersion interactions between two apposed graphene layers, a graphene 
layer and a substrate, and in a multilamellar graphene system are analyzed within the framework of 
the Lifshitz theory. This formulation hinges on a known form of the dielectric response function of 
an undoped or doped graphene sheet, assumed to be of a random phase approximation form. In the 
geometry of two apposed layers the separation dependence of the van der Waals-Casimir interaction 
for both types of graphene sheets is determined and compared with some well known limiting cases. 
In a multilamellar array the many-body effects are quantified and shown to increase the magnitude 
of the van der Waals-Casimir interactions. 
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I. INTRODUCTION 

Graphene appears to be the only known mono-atomic 
two-dimensional (2D) crystal and apart from the intrin- 
sic interest it engenders, it is becoming more and more 
also a focus of possible and desired advanced technologi- 
cal applications^. It is for these reasons that in the past 
several years we have witnessed a veritable explosion of 
theoretical and experimental interest in graphene^. The 
Nobel prize for physics in 2010 only consolidated this 
trend. Graphene differs fundamentally from other known 
2D semiconductors because of its unique electronic band 
structure, viz. the monoatomic sheet of carbon atoms ar- 
ranged in a honeycomb lattice leads to an electron band 
structure that displays quite unusual properties^. The 
Fermi surface is reduced to just two points in the Bril- 
louin zone and the value of the band gap is reduced to 
zero. The energy dispersion relation for both the con- 
duction and the valence bands are linear at low energy, 
namely less than 1 eV, meaning that the charge carri- 
ers behave as relativistic particles with zero rest mass. 
The agent responsible for many of the interesting elec- 
tronic properties of graphene sheets is the non-Bravais 
honeycomb-lattice arrangement of carbon atoms, which 
leads to a gapless semiconductor with valence and con- 
duction TT-bands. 

States near the Fermi energy of a graphene sheet are 
described by a massless Dirac equation which has chi- 
ral band states in which the honeycomb-sublattice pseu- 
dospin is aligned either parallel or opposite to the en- 
velope function momentum. The Dirac-like wave equa- 
tion leads to both unusual electron-electron interaction 
effects and to unusual response to external potentials. 
When the graphene sheet is chemically doped with ei- 
ther acceptor or donor impurities its carrier mobility can 
be drastically decreased^. Because of its 2D periodic 
structure graphene is closely related to single wall car- 
bon nanotubes, being in fact a carbon nanotube rolled 



out into a single 2D sheet^. The main difference between 
the electronic properties of single wall carbon nanotubes 
and graphene is that the former show circumferential pe- 
riodicity and curvature that leave their imprint also in 
the electronic spectrum and consequently also van der 
Waals (vdW) interactions^. 

On the other hand, graphite appears to be the poor 
cousin of graphene though it is the stable form of car- 
bon at ordinary temperatures and pressures. Many ef- 
forts have been invested into understanding its structural 
and electronic details (for an account see Ref4) . Various 
known modifications of graphite differ primarily in the 
way the mono-atomic two-dimensional graphene layers 
stack. Their stacking sequence in terms of commonality 
is ABA for the Bcrnal structure, AAA for simple hexag- 
onal graphite or ABC for the rhombohedral graphite^,. 

Graphene layers in graphitic systems are basically 
closed shell systems and thus have no covalent bond- 
ing between layers which makes them almost a perfect 
candidate to study long(er) ranged non-bonding interac- 
tions. Indeed, they are stacked at an equilibrium inter- 
layer spacing of about 0.335 nm and are held together 
primarily by the non-bonding long range vdW interac- 
tionsiS. Therefore the interaction between graphene lay- 
ers can be described as a balance between attractive vdW 
dispersion forces and corrugated repulsive (Pauli) overlap 
forcesii, following in this respect closely the paradigm of 
nano-scale interactional^. 

Besides a few notable exceptions^, until 2009 many 
electronic and optical properties of graphene could be ex- 
plained within a single-particle picture in which electron- 
electron interactions are completely neglected. The dis- 
covery of the fractional quantum Hall effect in graphenei^ 
represents an important hallmark in this context. By now 
there is a large body of experimental worki^"— showing 
the relevance of electron-electron interactions in a num- 
ber of key properties of graphene samples of sufficiently 
high quality. 
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Because of band chirality, the role of electron-electron 
interactions in graphene sheets differs in some essential 
waysi^"— from the role which it plays in an ordinary 2D 
electron gas. One important difference is that the contri- 
bution of exchange and correlation to the chemical poten- 
tial is an increasing rather than a decreasing function of 
carrier density. This property implies that exchange and 
correlation increases the effectiveness of screening, in con- 
trast to the usual case in which exchange and correlation 
weakens screening^!. This unusual property follows from 
the difference in sublattice pseudospin chirality between 
the Dirac model's negative energy valence band states 
and its conduction band statesi^ii^, and in a uniform 
graphene system is readily accounted for by many-body 
perturbation theory. 

In this work we focus our efforts on the vdW dispersion 
component of the graphene stacking interaction. Dis- 
persion forces can be formulated on various levels^ giv- 
ing mostly consistent results for their strength and sep- 
aration dependence. In the context of graphene stack- 
ing interactions, the problem can be decomposed into 
the calculation of the dielectric response of the carbon 
sheets and the subsequent calculation of the vdW interac- 
tions either via the quantum-field-theory-based Lifshitz 
approach, as advocated in this paper, by means of the 
electron correlation energy2^ or the non-local vdW den- 
sity functional theory^l. One can show straightforwardly 
that in fact the non-local van der Waals functional ap- 
proach of the density functional theory and the Lifshitz 
formalism are in general equivalenti^^. 

Specifically we will calculate the vdW-Casimir inter- 
action free energy, per unit area between two graphene 
sheets as a function of the seperation between them, in 
a system composed of 

a) two apposed undoped or doped graphene sheets, 

b) an undoped or a doped graphene layer over a semi- 
infinite substrate, and 
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FIG. 1: (Color online) Schematic presentation of two 
graphene sheets of finite thickness immersed in vacuo at a 
separation of h. The thickness of both (left and right) layers 
are equal to a. In view of later generalizations, we have la- 
beled the left semi-infinite vacuum layer with L, the graphene 
layers with A, the intervening vacuum layer with B and the 
right vacuum layer with R. The form of the dielectric func- 
tions of the graphene layers is given in Eq. (|8]). 



c) a multilayer (infinite) array of graphene sheets. 

In the latter case we will investigate the many-body non- 
pairwise additive effects in the effective interaction be- 
tween two sheets within a multilayer array. We should 
note that non-pairwise additive effects are ubiquitous in 
the context of vdW interactions^^ often leading to non- 
trivial properties of macromolecular interactions. In this 
case they will lead to variations in the equilibrium stack- 
ing separation as a function of the number of layers 
in a graphitic configuration. In the calculation of the 
vdW-Casimir free energy we will employ the dielectric 
response function of a single graphene layer calculated 
previouslyi^. 



II. VDW-CASIMIR INTERACTION BETWEEN 
TWO LAYERS OF GRAPHENE 

The geometry of the system composed of two parallel 
graphene layers with thicknesses a, facing each other in a 
bilayer arrangement at a separation 6, is shown schemat- 
ically in Fig. [TJ In view of later generalizations we la- 
bel the left semi-infinite vacuum space as L, graphene 
sheets as A, the intervening layer as B and the right 
semi-infinite vacuum space by R. 

In order to calculate the vdW-Casimir dispersion in- 
teraction free energy in the planar geometry we use the 
approach of Refi^ where it has been calculated exactly 
for a multilayer planar geometry. The thus derived gen- 
eral form of the interaction free energy per unit area in- 
cludes retardation effects and is therefore valid for any 
spacing between the layers. 

For the system which is shown schematically in Fig. [TJ 
the vdW-Casimir interaction free energy is obtained in 
the Lifshitz form as 
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Q 



n=0 



(1) 



where Fgg{b) stands for the graphene-graphene interac- 
tion free energy as a function of the layer spacing b (nor- 
malized in such a way that it tends to zero at infinite 
inter layer separation). In the above Lifshitz formula the 
Q summation is over the transverse wave vector and the 
n summation (where the prime indicates that the n = 
term has a weight of 1/2) is over the imaginary Matsub- 
ara frequencies 
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where ks is the Boltzman constant, T is the absolute 
temperature, and h is the Planck constant devided by 
27r. All the quantities in the bracket depend on Q as 
well as 

The other quantities entering the Lifshitz formula are 
defined as 
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with 
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where Ai i_i quantifies the dielectric discontinuity be- 
tween homogeneous dielectric layers in the system, where 
a layer labeled by i — 1 is located to the left hand side of 
the layer labeled by i (for details see Refi^). Also Ki{t£,n) 
for each electromagnetic field mode within the material 
i is given by 



(5) 



where c is the speed of light in vacuo, Q is the magnitude 
of the transverse wave vector, and ei{t^n) and /ii(i^„) 
are the dielectric function and the magnetic permeability 
of the i-th layer at imaginary frequencies, respectively. 
For the sake of simplicity we assume that for all layers 
/^j(*^n) = 1 and that the dielectric function for vacuum 
layers equals to 1 for all frequencies. 

Note that ei(i^) is standardly referred to as the vdW- 
London transform of the dielectric function and is defined 
via the Kramers-Kronig relations as^i 
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It characterizes the magnitude of spontaneous electro- 
magnetic fluctuations at frequency ^. In general e(«^) is 
a real, monotonically decaying function of the imaginary 
argument ^ (for details see Parsegian's book in Ref.— ). 

In order to proceed one needs the vdW-London trans- 
form of the dielectric function of all layers in the system. 
The detailed Q- and w-dependent form of the dielectric 
function for undoped and/or doped graphene layers are 
introduced in Sees. Ill A\ and HTbI 



A. Two undoped graphene layers 

We employ the response function of a graphene layer 
from Refs i^^i^^ which for doped graphene assumes the 



form 

x(Q,<«,Ai 7^ 0) 



=Re 



^ 2/i + z^„;t ^ 
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vQh 
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where g = A, v ^ 10^ m/s is the Fermi velocity in 
graphene layer and ^ = ep = hvkp is the chemical po- 
tential, ep the Fermi energy, and hkp is the Fermi mo- 
mentum, where kp = [Airp/ g)^/"^ and p is the the average 
electron density. 

To begin with, we assume that two layers are decoupled 
and ignore the interlayer Coulomb interaction. The vdW- 
London dispersion transform of the dielectric function on 
the level of the random phase approximation (RPA) is 
then given by 



e(Q,ze„) = l-y(Q)x(Q,«^n,//7^0), 



(8) 



where V{Q) is the (transverse) 2D Fourier-Bessel trans- 
form of the Coulomb potential, V{Q) ~ 4^eoe Q ■ ^ ^® 
electric charge of electron, eo is the permittivity of the 
vacuum and is the average of the dielectric constant 
for the surrounding media which is equal to 1 for vacuum. 
In what follows we furthermore assume that to the low- 
est order the dielectric properties of the graphene layers 
are not affected by the variation of the separation be- 
tween them. This assumption is also consistent with the 
Lifshitz theory that presumes complete independence of 
the dielectric response functions of the interacting layers. 

For undoped graphene layer, i.e. p = 0, the expres- 
sion for the vdW-London dispersion transform of the di- 
electric function simplifies substantially and assumes the 
form 



e((3,<r,) = 1 + 



iragcQ 



(9) 



where a is the electromagnetic fine-structure constant 

2 1 

a = . ^ > ~ It should be noted that, in general, 

47reo'ic 137 ' ^ ' 

a model going beyond the RPA is necessary in order to 
account for enhanced correlation effects that would be 
present in an undoped system^^. In this paper, however, 
we restrict ourselves to the RPA approximation and an- 
alyze its predictions in detail. 
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FIG. 2: (Color online) Magnitude of the interaction free en- 
ergy per unit area of the system composed of two undoped 
graphene layers immersed in vacuo at an interlayer spacing 
6, and at temperature 300 K. The functional dependence 
of the interaction free energy on 6 is compared with the fol- 
lowing scaling forms: and as the separation 
increases. 



The functional dependence of the interaction free en- 
ergy of the system per unit area, Eq. ([T}, is presented 
in Fig. [5] as a function of the separation between two 
graphene layers. We assume that the graphene layers 
are immersed in vacuo and both of them have the same 
thickness lA as well as equal susceptibilities. Note that 
in all cases considered in this paper the interaction free 
energies as defined in Eq. ([1]) are negative reflecting at- 
tractive vdW-Casimir force between graphene layers in 
vacuum. For the sake presentation, we shall plot the ab- 
solute value (magnitude) of the free energy in all cases. 

As one can discern from Fig. [5] the general dependence 
of the vdW-Casimir interaction free energy on the sepa- 
ration between the graphene layers has the scaling form 
of a power law, 5^", with a weakly varying separation- 
dependent scaling exponent, n(6). This scaling exponent 
can be defined standardly as^ 



n[h) 



dlnFggjb) 
dlnb ' 



(10) 



For two undoped graphene sheets we observe that at 
small separations the functional dependence of the free 
energy on interlayer spacing yields the scaling exponent 
n = 3 for smallest values of the separation. The scaling 
exponent then steadily increases to n = 4, then n ~ 5, 
and finally at asymptotically large separations it reverts 
back to n = 4. This variation in the scaling exponent for 
the separation dependence of the interaction free energy 
can be rationalized by invoking some well known results 



on the vdW interaction in multilayer geometries (see e.g. 
the relevant discussions in Refi^). 

For example, for two semi-infinite layers the interac- 
tion free energy should go from the non-retarded form 
characterized by n = 2 for small spacings, through re- 
tarded n — S form for larger spacings and then back to 
zero-frequency-only form that also scales with n = 2 but 
with a different prefactor than the non-retarded form. 
For two infinitely thin sheets, on the other hand, we 
have the non-retarded n = 4 form for small separation, 
followed by the retarded n = 5 form for larger spacings 
and then reverting back to zero- frequency-only term with 
n = A scaling, but again with a different prefactor than 
the non-retarded limit. Furthermore, the transitions be- 
tween various scaling forms and the locations of the tran- 
sition regions are not universal but depend crucially on 
the characteristics of the dielectric spectra and can thus 
be quite complicated, sometimes not yielding any easily 
discernible regimes with a quasi-constant scaling expo- 
nent n. 

Reading Fig. [5] with this in mind we can come up with 
the following interpretation of the calculated separation 
dependence: for small separation the interaction free en- 
ergy is dominated by the n = 4 dependence except in the 
narrow interval very close to vanishing separation where 
the dependence asymptotically levels off at n < 4 form. 
This form is consistent with the non-retarded interaction 
between two very thin layers, see above, for small but 
not vanishing separations. For vanishing interlayer sep- 
arations the final leveling-off of the scaling exponent is 
due to the fact that the system is approaching the limit 
of two semi- infinite layers where in principle n — > 2, but 
in reality the finite thickness of the graphene sheets is 
way too small to observe this scaling in its pure form. 
All we can claim is that for vanishing interlayer spacings 
the scaling exponent drops below the n ~ A value, valid 
for two infinitely thin layers. 

For larger values of the interlayer separations we then 
enter the retarded regime with n = 5 scaling exponent, 
again valid strictly for two infinitely thin layers. The re- 
tarded regime finally gives way to the regime of asymp- 
totically large spacings where the interaction free energy 
limits towards its form given by the zero frequency term 
in the Matsubara summation and characterized by n = 4 
scaling dependence. Obviously the numerical coefficient 
in the small separation non-retarded and asymptotically 
large separation regimes (both with n = 4) are necessar- 
ily different. 

The interaction free energy scaling with the interlayer 
separation is thus completely consistent with the vdW- 
Casimir interactions between two thin dielectric layers for 
all, except for vanishingly small, separations where finite 
thickness effects of the graphene sheets leave their mark 
in a smaller value of the scaling exponent that should ide- 
ally approach the value valid for a regime of interaction 
between two semi-infinite layers. 

The numerical value of the interaction free energy per 
unit area at 1 nm is about 5.64 x 10^^ eV/m^. That 
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means that for two graphene layers with surface area 
of 10~^^ the magnitude of the free energy is about 
564 eV at 1 nm separation; at 10 nm it is about 0.01 eV 
and at 100 nm about 3.6 X lO^'^ eV for the same surface 
area. 



B. Two doped graphene layers 

For a doped graphene layer the vdW-London disper- 
sion transform of the dielectric function can be read off 
from Eqs. ^ and Q as 

e(Q,<«,/iy^O) = 1 
nagcQ 



, i arcsm [-Ai - -Bi\ 

+2^1^ {Aj + Bl) cos [\ arg(A2 + ^S^)] 

-^{Al + Bl)"\in [\ arg(^ + }, (11) 
with the following coefficients 



Ai 



vQ 



^2 = 1+1^1 - 16 

B2 = 



np 



vQj "gQ^' 

£.n\/4.TTp/g 



(12) 



With this dielectric response function we again evalu- 
ate the vdW-Casimir interaction free energy of the sys- 
tem per unit area, Eq. ([1]), as a function of the separation 
between two graphene layers b as shown in Fig. [31 The 
interaction free energy has again a scaling form with a 
scaling exponent varying with the separation between the 
layers. It is clear that in this case it is much more dif- 
ficult to partition the variation of the scaling exponent 
into clear-cut piecewise constant regions. 

As it can be seen from Fig. [3] at small separations the 
form of the functional dependence of the interaction free 
energy has n < 3. Then for increasing spacings there 
follows a relatively broad regime with n = 3 — 4, fol- 
lowed eventually by the scaling form with n = 2 for 
6 > 5 X 10~®m. One needs to add here that only the 
scaling regime of n = 2 for asymptotically large sepa- 
rations and an intermediate regime with n = 3 — 4 arc 
clearly discernible. 
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FIG. 3: (Color online) Magnitude of the interaction free 
energy per unit area of the system composed of two doped 
graphene layers (solid curve) compared with that of two ideal 
metallic sheets (dashed curve) immersed in vacuo as a func- 
tion of the interlayer spacing 6, at temperature 300 K. The 
functional dependence of the free energy on the interlayer 
spacing is compared with scaling forms b~^, and 

in various regimes of separation. 



We can gain some understanding of these regimes by 
comparing with the various exact limits in the layer ge- 
ometry as before. Such comparison is however not as 
straightforward as before. The asymptotic n = 2 regime 
is easiest to rationalize: it has the same scaling form as 
the finite temperature vdW-Casimir interaction between 
two metallic sheets at asymptotically large separations. 
The presence of free charges would in fact be a rea- 
sonable characterization of doped graphene layers. For 
smaller separations we then enter the regime dominated 
by the retardation effects with n ~ 3 — 4 form, and fi- 
nally for vanishing separations we approach the regime of 
n = I . Recent calculations of vdW interactions between 
thin metallic layers indeed lead to exactly this exponent 
for small layer separations—. The doped graphene sheet 
results would thus indicate that the dependence of the 
vdW-Casimir interactions free energy on the separation 
could be rationalized in terms of interactions between two 
thin metallic sheets. 



For comparison we have also plotted the interaction 
free energy between two ideal metallic sheets which ex- 
hibits a much stronger attractive interaction free energy, 
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FIG. 4: Magnitude of the free energy per unit area of the sys- 
tem composed of two undoped (dashed line) and doped (solid 
line) graphene layers (with the electron density p = 10^^ /m^) 
immersed in vacuo as a function of the interlayer spacing b, 
at temperature 300 K. As seen for all separations the magni- 
tude of the interaction free energy for doped graphene layers 
is greater than that of the undoped one. 



FIG. 5: (Color online) The rescaled interaction free energy, 
^aai'n) / Pgaiv ~ 0)i f™' the system composed of two doped 
graphene layers with different electron densities immersed in 
vacuum as a function of = p\+p2 different interlayer 
separations b — 1 nm, 10 nm, 10 /im, and 100 fim (from top 
to bottom) at temperature 300 K. Note that Fgg{rj = 0) has 
been calculated for pi — p2 ~ 10^^/m^. 
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The numerical value of the energy per unit area at 1 nm 
is about 1.77 X 10^^ eV/m^ which is equal to 1.77 x lO'' eV 
for surface area of 10^^^ m^; at 10 nm it is about 16.9 eV 
and at 100 nm it is about 0.01 eV for the same surface 
area. 



C. Doped vs. undoped graphene 

It is instructive to compare the interaction between two 
graphene layers in the undoped and doped cases. For this 
purpose we have plotted the interaction free energy of the 
system for both cases in Fig. H) The electron density in 
doped graphene is assumed to be p = lO^^/m^ (solid 
curve). The dashed curve is the interaction free energy 
of the system composed of two undoped graphene layers. 
As it can be seen, for all separations the magnitude of 
the interaction free energy for doped graphene layers is 
more than that of the undoped one (note again that the 



interaction free energies as defined in Eq. ([1} are nega- 
tive in the present case due to attractive vdW-Casimir 
force between graphene layers in vacuum) . At separation 
1 nm the vdW-Casimir interactions for doped graphene 
is about 30 times the magnitude of the interaction for 
undoped graphene, while at the separation of 10 nm this 
ratio is about 1600. This means that the attractive in- 
teraction between graphene layers is enhanced when the 
contribution of the electron density in the dielectric func- 
tion of the graphene layers is taken into account. This 
same trend was observed also in the work of Serneliua^^ 
and is clearly a consequence of the fact that the largest 
value of vdW-Casimir interactions is obtained for ideally 
polarizable, i.e. metallic layers. The closer the system is 
to this idealized case, the larger the corresponding vdW- 
Casimir interaction will be. 

Let us investigate also the effect of asymmetry of 
doped graphene sheets on vdW-Casimir interactions be- 
tween them. Introducing the dimensionless parameter 



P1-P2 
P1+P2 ■ 



where pi is the electron density of the 



■n 

th graphene layer (z = 1,2), we find that the interac- 
tion free energy depends on the asymmetry in the sys- 
tem. When 77 = 0, the electron densities arc the same 
for both graphene layers and we have a symmetric case, 
whereas 77 = 1 means that one of the graphene layers 
is undoped while the other one is doped, leading to an 
asymmetric case. In Fig. [S] we have plotted the rescaled 
free energy Fggiv)/ Pggiv — 0) of the system composed 
of two graphene layers as a function of 77 for different val- 
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FIG. 6: Hamaker coefficient for a system of two undoped 
graphene layers as a function of the distance between them, 
b. The top sohd hne shows the general form of the Hamaker 
coefficient as defined via Eq. p4|l . the dashed line shows the 
limiting large-distance form, Eq. ([ISp . and the dotted line 
shows the limiting small-distance form, Eq. 1)16^ . 



ues of the interlayer separation. The magnitude of the 
electron density for one of the layers has been fixed at 
pi = while that of the other layer, p2, varies. 

As seen from this figure, the curves show a monotonic 
dependence on 77 with a stronger interaction at smaller 
values of 77. Note that at large separations the curves 
tend to coincide and will become indistinguishable. The 
asymmetry effects are therefore largest at small separa- 
tions between the interacting graphene layers. 



D. Hamaker coefficient for two grapliene layers 

The general form of the Hamaker coefficient, Agg, for 
a system composed of two dielectric layers of finite thick- 
ness (Fig. [1]) when retardation effects are neglected is 
defined via3^ 
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At large separations b ^ a, the Hamaker coefficient, 
ained fi 



A^gg^'^. can be obtained from 
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FIG. 7: (Color online) Schematic presentation of a graphene 
layer of thickness a (labeled by A) apposed to a substrate 
(labeled by L) at a separation b. We have labeled the inter- 
vening layer (assumed to be vacuum) with B, and the right 
one with R. The dielectric function of the graphene layer is 
defined via Eq. ((SJ and for substrate via Eq. (|f 81) . 



In Fig. [51 we show the Hamaker coefficient as a func- 
tion of the layer separation for a system of two un- 
doped graphene layers. We show the general form of 
the Hamaker coefficient from Eq. (fT4)) (solid line) as well 
as the limiting forms at small (dotted line, Eq. ((T5)) ) 
and large (dashed line, Eq. (US])) separations. The 
large-distance limiting form obviously coincides with the 
general form at separations beyond 50 nm. The small- 
distance limiting form tends to the general form at small 
separations but given that the thickness of the layers is 
only about lA, it is expected to merge with the general 
form in sub- Angstrom separations. For doped graphene, 
a similar analysis to Eq. is not possible because the 
corresponding general expression for the interaction free 
energy is missing. 



III. VDW-CASIMIR INTERACTION BETWEEN 
A GRAPHENE LAYER AND A SEMI-INFINITE 
SUBSTRATE 

In this section we study the interaction between a 
graphene layer and a semi-infinite dielectric substrate as 
depicted schematically in Fig. [T] For this system the free 
energy per unit area is 



Fsgjb) 

s 



= knT X 



EE In 

Q n=0 



-2aK / 



1 + A^sAij^e-^'"^^ 



AsLe 



(17) 



(16) 



where Fag now stands for the interaction free energy be- 
tween the substrate and the graphene layer. We have 
excluded the explicit dependence of the quantities in the 
bracket on the imaginary Matsubara frequencies, but 
they are the same as in Eq. (|4]). 

In order to gain insight into the magnitude of the vdW- 
Casimir interaction free energy and for the sake of sim- 
plicity, we assume that the semi-infinite substrate is made 
of Si02 which has the vdW-London dispersion transform 
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FIG. 8: (Color online) Magnitude of the interaction free 
energy of the system composed of a Si02 substrate and an 
undoped graphene layer, Eq. (|17|) . plotted at temperature 
300 K as a function of the separation b. The functional depen- 
dence of the free energy on b is compared with scaling forms 
b~^, b~^, b~'^ and 6""^ in various regimes of separation. 



of the dielectric function of the for m-^^i^^ 



f-L{lin) = 1 + 



■'uv 



£2 



(18) 



'^IR 



where the values of the parameters, Cuv = 1-098, Cir = 
1.703, wuv = 2.033 x lO^*' rad/s, and ujyr, = 1.88 x 10^'' 
rad/s have been determined from a fit to optical data^. 
The static dielectric permittivity of Si02 is then obtained 
as e(0) ~ 3.81. A characteristic feature of the vdW- 
London transform of Si02 is thus that it contains two 
relaxation mechanisms. The first one is due to electronic 
polarization and the second one is due to ionic polariza- 
tion. All calculations of the vdW-Casimir interaction free 
energy are done at room temperature (300 K). 



A. Undoped graphene apposed to a substrate 

Using the vdW-London transforms of the dielectric 
functions given in the preceding sections, one can now 
calculate the free energy, Eq. (fT7|) . for an undoped 
graphene layer next to a semi-infinite Si02 substrate. 
The results are shown in Fig. [S] 

At small separations, the free energy varies with a scal- 
ing exponent n = 2, while at larger separations one can 
distinguish the scaling regimes n = 3, and ti = 4, fi- 
nally approaching the asymptotic limit at large separa- 
tions with n ~ 3. This sequence of interaction free energy 
scalings can be rationalized as follows: at asymptotically 
large separations we are at the zero- frequency Matsubara 



FIG. 9: (Color online) Magnitude of the interaction free en- 
ergy of the system composed of a Si02 substrate and a doped 
graphene layer, Eq. (|17|) is plotted at temperature 300 K as 
a function of the separation b. Here the functional depen- 
dence of the interaction free energy for the graphene doping 
electron density p = 10^® /m^ (solid line) is compared with 
scaling forms b~^ and b~^ in various regimes of separa- 

tion. We also include the same plot for the doping electron 
density p = 10^* /m^ (dashed line). 



term for a semi-infinite layer and a thin sheet. This case 
is right in between the asymptotically large separation 
limit for two semi-infinite layers (n = 2) and two in- 
finitely thin layers (n = 4). For smaller spacings we then 
progressively detect contributions from higher Matsubara 
terms which lead to scaling exponent n = A that corre- 
sponds to a retarded form of the interaction free energy 
and then for yet smaller spacings the scaling exponent re- 
verts to n = 3 non-retarded form of the vdW interaction 
between a semi-infinite substrate and an infinitely thin 
sheet. For vanishing spacings the finite thickness of the 
sheet starts playing a role and eventually, we approach 
the n = 2 scaling for two semi-infinite layers. 

The magnitude of the interaction free energy per unit 
area at 1 nm is about 1.57 x 10^^ eV/m^. It means that 
for a system with surface area of 10^^^ m^ this value is 
about 1.57 X 10^ eV. At 10 nm it is about 0.36 eV and 
at 100 nm it is about 8.71 x 10"^ eV for the same surface 



B. Doped graphene apposed to a substrate 

The vdW-Casimir interaction free energy Eq. ((TT)) for 
a system composed of a doped graphene layer next to 
a semi-infinite Si02 substrate is shown Fig. [5] for two 
values of the electron density p = lO^^'/m^ (dashed line) 
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FIG. 10: Hamaker coefficient for a system of an undoped 
graphene layer apposed to a Si02 substrate as a function 
of the distance between them, b. The top solid line shows 
the general form of the Hamaker coefficient as defined via 
Eq. (|19|) . the dashed line show the limiting large-distance 
form, Eq. (|20|) . and the dotted line shows the limiting small- 
distance form, Eq. (|2ip . 



and p = 10 °/m- (solid line). 

At small separations the free energy shows the n = § 
scaling, while for larger spacings it shows a scaling ex- 
ponent 71 = 3, approaching the n — >■ 2 limit for asymp- 
totically large separations. In this respect the case of a 
thin doped graphene sheet apposed to a semi-infinite sub- 
strate is very similar to the case of two thin doped layers, 
except that the retarded regime covers a smaller interval 
of spacings. This means that the metallic nature of one 
of the interacting surfaces is enough to switch the behav- 
ior of the interaction free energy completely towards the 
case of two metallic interacting surfaces. This case has 
in fact not yet been thoroughly discussed in the litera- 
ture. The changes in the slope appear to occur at the 
same values of the interlayer spacing when the electron 
density decreases (compare dashed and solid curves). 

The magnitude of the free energy for the electron 
density p = 10^"'/m^ (dashed line) and surface area of 
^2 ^-^ jg abotit 2.12 X 10^ eV, while at 10 nm 
it is about 3.28 eV, and at 100 nm is about 4.20 x IQ-^ cV. 

These values increase for larger electron densities, e.g., 
for p = 10^^/m^ (solid line) and the same surface area, 
the free energy magnitude is 5.77 x 10'^ eV at 1 nm, while 
it is about 21.1 eV at 10 nm and about 2.67 x lO'^ eV 
at 100 nm. 



C. Hamaker coefflcient for the graphene-substrate 
system 

The general form of the Hamaker coefficient, Agg, for 
a system composed of a dielectric layer of finite thickness 
apposed to a semi-infinite dielectric substrate (Fig. [7]) 
when retardation effects are neglected is defined via^ 



S 



1 - 



62 



{b + ay- 



(19) 



At large separations 6 ^ a, the Hamaker coefficient, 
lined fr 

Fsgib) 



A^^gS"^, can be obtained from 



(20) 



S 67r63 ' 

while at small separations 6 ^ a, it can be read off from 
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In Fig. [ini we again show the Hamaker coefficient as a 
function of the layer separation for a system comprising 
a Si02 substrate and an undoped graphene layer. We 
show the general form of the Hamaker coefficient from 
Eq. (fT9|) (solid line) as well as the limiting forms at large 
(dashed line, Eq. and small (dotted line, Eq. (PT|) ) 

separations. The large-distance limiting form obviously 
coincides with the general form at separations beyond 
10 nm but again since the thickness of the graphene layer 
is only about lA, the small-distance limiting form is ex- 
pected to merge with the general form in sub- Angstrom 
separations. 



IV. VDW-CASIMIR INTERACTION IN A 
SYSTEM COMPOSED OF N + 1 LAYERS OF 
GRAPHENE 

In this section we shall use the Lifshitz formalism in 
order to study the many-body vdW interactions in a sys- 
tem composed of -I- 1 layers of graphene. The layers 
are separated from each other by N layers of vacuum and 
are bounded at the two ends by two semi-infinite dielec- 
tric slabs as depicted in Fig. [TT] The thickness of each 
graphene layer is a while the separation between two suc- 
cessive layers is b. We have labeled the left semi-infinite 
dielectric medium with L, the right one with R (which 
both will be assumed to be vacuum) , the graphene layers 
with A and the vacuum layers with B. Following Refi^, 
one can calculate the vdW-Casimir part of the interaction 
free energy, -Fat (a, b), in an explicit form for any finite N. 
Interestingly, it turns out that for very large values of 
the vdW-Casimir free energy can be written as a linear 
function of N so that the interaction free energy FN{a, b) 
becomes^ 



FNia,b) = Nfgg{a,b) 



N:s>i, 



(22) 
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FIG. 11: (Color online) Schematic picture of the system com- 
posed of A'' + 1 layers of graphene (^4) with equal thicknesses 
a separated from each other by layers of vacuum (B) with 
equal thicknesses of b. The two semi-infinite substrates in the 
left and right end of the system are labeled by L and R and 
will be assumed to be vacuum. 
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where fgg(a,b) can be interpreted as an effective pair 
interaction between two neighboring layers in the stack 
and is given by 



fgg{a, h) = ksT X 



n=Q 



(1 - A2e-2«^") 



G{a,b,A) 



y (1 - A2e-2«^'i)2 
with Gia, b, A) defined as 



FIG. 12: (Color online) Magnitude of the interaction free 
energy per unit area and number of layers, \fgg\/S, plotted 
as a function of the separation, 6, between two successive un- 
doped (bottom black dot-dashed line) and doped (top black 
solid line) graphene layers for a system composed of infinitely 
many layers of undoped/doped layers as schematically de- 

picted in Fig. 1111 The temperature of the system is 300 K, 

the thickness of the graphene layers is fixed at lA and we 
have chosen the dielectric function for the undoped case from 
Eq. ((9| and for doped case (with the doping electron den- 
(23) sity chosen as p = 10^^/m^) from Eq. (fTTj) . The functional 



4e- 



-2{KAa+KBb) 



(24) 



Here A is 



where ka and kb are 



KA^B - I^B^A 
KA^B + I^B^A ' 



(25) 



(26) 



For simplicity we have dropped the explicit dependence 
on the Matsubara frequency in all the above expressions. 
We stress that fgg{a,b) is an effective pair interaction 
between two neighboring graphene layers in a multilayer 
geometry and is in general not equal to Fggia^b) in Eq. 
([T]) , which is valid for two interacting layers in the absence 
of any other neighboring layers. The difference between 
these two interaction free energies thus encodes the non- 
pairwise additive effects in the interaction between two 
layers due to the presence of other vicinal layers. 



form of the free energy for the undoped case is compared with 
scaling forms , and &~* in various regimes of sep- 

aration. For the doped case, the free energy is compared with 
the scaling forms h^^'^ and h~'^ . 



A. A'^ + 1 undoped graphene layers 

Let us first consider the case of iV+1 undoped graphene 
layers. In this case the vdW-Casimir interaction free en- 
ergy per unit area and per number of layers, \fgg{b)\/ S 
(black dot-dashed line), has been plotted in Fig. [T2l as a 
function of the separation between the layers, b. The tem- 
perature of the system is chosen as 300 K, the thickness 
of the graphene layers is 1 A and we have used the dielec- 
tric function given by Eq. ^ for each undoped graphene 
sheet. 

The value of \fgg\/S at 1 nm is about 5.9 x lO^^ eV/m^ 
which is 5.9 x 10^ eV when the surface area is equal to 
10~^2 m^. At 10 nm the value of \fgg\ for the same 
surface area is about 0.01 eV and at 100 nm is about 
3.9 X 10-^ eV. 

The scaling of \fgg{b)\/S for different values of the in- 
terlayer spacing b is shown in Fig.[T3 It shows the scaling 
exponent n = 3 at vanishing separations while at finite 
yet small separations it is characterised by ri = 4, contin- 
uously merging into a n = 5 form and finally attaining 
the n = 4 form. The rationalization of this sequence of 
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FIG. 13: (Color online) Magnitude of the interaction free 
energy per unit area and number of layers, \fgg\/S, for the 
system composed of A'^ + 1 layers of undoped (bottom black 
dot-dashed line) and doped (top black solid line) graphene is 
compared with the magnitude of the interaction free energy 
per unit area, \Fgg\/S, for the system composed of only two 
undoped (green dots) and doped (red dots) graphene layers 
as a function of the separation between the layers, b. The 
inset shows the ratio of these two quantities for both undoped 
(black dot-dashed line) and doped (black solid line) cases. 



scaling exponents is exactly the same as in the case of 
two isolated layers and will thus not be repeated here. 

One can directly compare the reduced free energy, 
fgg{b), with that of a system composed of only two un- 
doped graphene layers of the same thickness a, Fggih) 
(i.e., comparing the results in Fig. [2] with the correspond- 
ing results in Fig. [T2|) . This is shown in Fig. [13] (black 
dot-dashed line and green dots), where apparently the 
results nearly coincide. However, by inspecting the ratio 
between these two interaction free energies it turns out 
that in the multilayer system the interaction free energy 
per layer is slightly more attractive than in the case of 
a two-layer system. This difference thus stems directly 
from the many-body effects which in this case augment 
the binding interaction in a graphitic stack. 



B. + 1 doped graphene layers 

In Fig. [T2]the magnitude of the vdW-Casimir interac- 
tion free energy for doped graphene layers per unit area 
and per number of layers, \fgg{b)\/S, is plotted (black 
solid line) as a function of the separation between lay- 
ers, b. The vdW-London dispersion transform of the 
dielectric function for each graphene sheet is chosen as 
in Eq. pT|) . We have fixed the density of electrons for 



all the graphene layers as p = 10^^/m^. The value of 
\fgg\/S at 1 nm is about 1.96 x 10^^ eV/m^, which is 
about 1.96 X 10^ eV when the surface area is equal to 
10~i2 m^. At 10 nm the value of \fgg\ is about 18.7 eV 
and at 100 nm it is about 1.38 x 10^^ eV. The scaling 
exponents of the interaction free energy dependence for 
different regions of interlayer spacings are illustrated in 
Fig. [T^l The scaling exponent is n = 2.5 for small sep- 
arations while at larger separations it tends towards the 
value n = 3. It is thus exactly the same as in the case 
of two isolated doped graphene sheets, see Fig. [31 except 
that in the multilayer geometry we have not shown the 
same range of separations as for two isolated layers. 

The comparison of vdW-Casimir interaction free ener- 
gies in the case of two isolated doped graphene sheets 
with the effective interaction between two graphene 
sheets in a multilayer system (i.e., comparing \Fgg\ in 
Fig. [3|with \ fgg \ in Fig. [12]) is made in Fig. [13] As shown 
by the inset the interaction free energy is again slightly 
more attractive within a multilayer. Comparing the re- 
sults in the inset of Fig. [13] shows that in average the 
many-body effects are stronger in the doped multilayer 
case than in the undoped case. 

Note also that a direct comparison between the un- 
doped and doped systems (Fig. [T^]) shows that for all 
separations the free energy magnitude for a doped mul- 
tilayer (solid line) is more than that of the undoped one 
(dot-dashed line) and thus the interaction is more attrac- 
tive in the former case. At separation 1 nm the magni- 
tude of the doped free energy is about 34 times larger, 
while at the separation 10 nm it is about 1730 larger than 
that of the undoped one. 



V. CONCLUSION 

In this work we have studied the vdW-Casimir inter- 
action between graphene sheets and between a graphene 
sheet and a substrate. We calculated the interaction free 
energy via the Lifshitz theory of vdW interactions that 
takes as an input the dielectric functions, or better their 
vdW-London transform, of isolated layers. Within this 
approach it would be inconsistent to take into account 
any separation dependent coupling between the dielectric 
response of the layers. This need possibly not be the case 
for some other approximate approaches to vdW interac- 
tions as in, e.g., the vdW augmented density functional 
theory (see the paper by Langreth ct al. in Ref<^). 

By inserting the random phase approximation dielec- 
tric function of a graphene layer into the Lifshitz theory 
we are thus in a position to evaluate not only the pair in- 
teraction between two isolated graphene sheets, but also 
between a graphene sheet and a semi-infinite substrate 
of a different dielectric nature (Si02 in our case) as well 
as the effective interactions between two graphene sheets 
in an infinite stack of graphene layers. All these cases 
that have been analyzed and discussed above are relevant 
for many realistic geometries in nano-scale systemsi^ and 
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thus deserve to be studied in detail. 

In the three cases studied we found the following salient 
features of the vdW-Casimir interaction dependence on 
the separation between the interacting bodies: 

1- In a system composed of two graphcne layers we 
demonstrated that the vdW-Casimir interactions 
in the case of undoped graphene show scaling ex- 
ponents identical to those displayed in the case 
of interacting thin dielectric layers. In the doped 
case the scaling exponents are consistent with vdW- 
Casimir interactions between two thin metallic lay- 
ers. 

2- In a system composed of a semi-infinite dielectric 
substrate and an undoped graphene layer the vdW- 
Casimir interactions display scaling exponents ex- 
pected for this asymmetric geometry. For a doped 
graphene layer the exponents revert to the previous 
case of two doped graphene layers. 

3- In a multilayer system composed of many graphene 
sheets the vdW-Casimir interaction scaling expo- 
nents are the same as in the case of two isolated 
layers but the interactions are stronger due to many 
body effects as a consequence of the presence of 
other layers in a stack. 

In order to describe the correlation effects especially 



at low doping or the interlayer coupling on a more sys- 
tematic level, one needs to go beyond the standard ran- 
dom phase approximation by incorporating more so- 
phisticated theoretical models for the dielectric response 
function which would be worth exploring further in the 
future. 

The main motivation for a detailed study of vdW- 
Casimir interaction between graphene sheets in graphite- 
like geometries is the fact that graphitic systems be- 
long to closed shell systems and thus display no covalent 
bonding, so that any bonding interaction is by neces- 
sity of a vdW-Casimir type. Its detailed characterization 
is thus particularly relevant for this quintessential nano- 
scale systemic. 
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